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Abstract 

In this paper we explore the large N limit of the glueball mass spectrum for 2+1 dimensional 
pure gauge theory. We employ Hamiltonian lattice gauge theory (LGT) and analytic variational 
techniques to calculate glueball masses for finite values of N < 25. The results are then accurately 
extrapolated to infinite N. An empirical observation is discussed which links the infinite N glueball 
spectrum to a simple oscillator model. 
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I. INTRODUCTION 

In this paper we study the glueball mass spectrum for large iV pure gauge theory (no 
quarks) in 2+1 dimensions using Hamiltonian LGT. Our primary motivation lies in recent 
developments in string theory where a direct connection between pure gauge theory in the 
large N limit and certain string theories has been conjectured. An independent test of 
this conjecture is provided by large iV calculations of the glueball mass spectrum in LGT. 
Impressive Monte Carlo calculations by Teper [1] and Lucini and Teper [2, 3] have progressed 
to the stage where accurate values of glueball masses have been calculated up to iV = 6 in 
2+1 dimensions. This has allowed an extrapolation to the N — > oo limit. With the analytic 
techniques of the Hamiltonian approach [4] it is possible to extend these calculations to 
much larger N, allowing a more reliable N — > oo extrapolation. 

The outline of this paper is as follows. We start in Section II with a historical background, 
where the original motivations for large N physics are discussed and the recent developments 
in string theory introduced. In Section III we use the analytic variational approach discussed 
in Ref. 4 to calculate glueball masses at finite values of N up to N = 25 in some cases. These 
results are then extrapolated to explore the large N limit of the 2+1 dimensional glueball 
mass spectrum. The N — > oo limits obtained here are compared to the results of the 
analogous Monte Carlo studies in the Lagrangian approach. 

II. BACKGROUND 

The nonperturbative region of quantum chromodynamics (QCD) has been examined 
predominantly with numerical techniques. While the numerical Monte Carlo simulations 
of Lagrangian LGT have made significant progress, there have been few developments in 
analytic techniques. One analytic technique which has received considerable attention in 
recent years is the large N limit. 

In 1974 't Hooft proposed a study of SU(iV) gauge theories in the large N limit rather 
than the physically interesting N = 3 case [5]. It was hoped that the SU(A) theory could be 
solved analytically and would be, in some sense, close to SU(3). Indeed it was shown that 
the large N theory simplifies drastically [5] but still captures at least some of the complexity 
of the SU(3) theory [6, 7]. Motivated by the fact that the QCD coupling constant g 2 is 
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a poor expansion parameter because its value depends on the energy scale of the process 
under consideration, t' Hooft considered an expansion in another dimensionless parameter 
of SU(iV) QCD, l/N. Based on an ingenious organisation of Feynman diagrams, 't Hooft 
was able to show that in the large N limit, keeping the 't Hooft coupling g 2 N fixed, only 
planar diagrams remain. While the remaining planar theory can be solved exactly in two 
dimensions, the three and four dimensional theories have not yet been solved. 

Recent developments in string theories have attempted to address this problem. In 1997 
Maldacena conjectured a correspondence, in the duality sense, between superconformal field 
theories and string theory propagating in a non-trivial geometry [8]. The idea behind du- 
ality is that a single theory may have two (or more) descriptions with the property that 
when one is strongly coupled the other is weakly coupled. A technique for breaking the 
conformal invariance and supersymmetry constraints of Maldacena's original proposal was 
later developed by Witten [9]. This led to the hope that nonperturbative pure SU(iV) gauge 
theories could be described analytically in 3+1 dimensions by their string theory dual. To 
develop this hope into a concrete solution would require perturbative expansions within the 
dual string theory. Such developments would seem to be some way off. 

In order to test the viability of using string theories to probe nonperturbative QCD 
independent tests are needed. The glueball mass spectrum of QCD provides a perfect 
laboratory. Recent Monte Carlo calculations have provided stable estimates of glueball 
masses up to N = 6 [1, 3] in 2+1 dimensions and up to iV = 5 [3] in 3+1 dimensions. From 
these estimates an iV — > oo limit can be extrapolated in each case. Convincing evidence 
of 0(1/N 2 ) finite N corrections have been demonstrated in 2+1 dimensions, verifying a 
specific prediction of 't Hooft's l/N expansion in the quarkless case. From the string theory 
side, present calculations require a strong coupling limit (N, g 2 N — > oo) to be taken. In 
this limit the string theory reduces to classical supergravity and the calculation of glueball 
masses, m(J PC ), is straightforward. In this limit the glueball is represented by a dilaton 
field whose mass can be extracted by solving the dilaton wave equation [10-12]. Despite the 
extreme approximations required surprising quantitative agreement with the weak coupling 
results from LGT were obtained. Later it was realised that, in 2+1 dimensions, lower mass 
states can be constructed from the graviton field [13]. With this approach the quantitative 



3 



agreement with LGT is lost but qualitative agreement in the form of the prediction 



m(0 ++ ) < m(2 ++ ) < m(r + ) 



(1) 



still holds. The qualitative agreement between string theories and LGT persists in 3+1 
dimensions [14] with 



for both theories. It appears that the identification of glueball states remains a problem. A 
serious difficulty is the removal of unwanted states which have masses of the same magnitude 
as the glueballs. Recipes for removing these unwanted states have been proposed [15, 16] 
but new parameters need to be introduced. 

The exploration of large N gauge theories outside of string theory has not been restricted 
to traditional Monte Carlo methods. Recent studies using the light-front Hamiltonian of 
transverse LGT have shown agreement with the Monte Carlo calculations of Teper [17]. In 
this approach explicit calculations at iV — > oo are possible without the need for extrapola- 
tion. Continuum calculations, without the use of string theory, have also commenced. An 
impressive series of papers [18-20] has led to specific predictions for the string tension for 
all N in 2+1 dimensions [20]. These predictions lie within 3% of lattice calculations up to 
N = 6 beyond which LGT results are not available. A recent review is available in Ref. 21. 

III. EXTRAPOLATION 
A. Introduction 

It this section we calculate 2+1 dimensional SU(A^) glueball masses for various iV using 
the analytic variational technique of Ref. 4. Our intention is to explore the large N limit 
of glueball mass spectrum by extrapolating the results obtained to their N oo limit. 
Employing analytic methods we can extend the calculation of the lowest energy massgaps to 
SU(25) on a desktop computer. This should be compared to current results in the Lagrangian 
approach which are limited, at present, to SU(6). With large values of N available a more 
reliable extrapolation to the N — > oo limit is possible. 

To simulate the ground, or (perturbed) vacuum, state with energy E , we use the one 



m(0 ++ ) < m(2 ++ ) < m(Q- + ), 



(2) 
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plaquette trial state, 



\<k) = exp 



(c £ ReQ) |0) 

\ X,l<j / 



(3) 



Here, |0) is the strong coupling vacuum defined by £"(x)\0) = for all %, x and a = 
1,2, ... , N 2 — 1. Sf{x) is the lattice chromoelectric field on the directed link running from 
the lattice site labelled by x to that labelled by x + at. Here i is a unit vector in the i 
direction. The directed square denotes the traced ordered product of link operators, Ui(x), 
around an elementary square, or plaquette, of the lattice, 



where a is the lattice spacing. 

As discussed in Ref. 4, the precise form of c((3), where (3 = N/g 2 , is unimportant in a 
scaling region. As it turns out, to calculate variational wave functions, by minimising the 
vacuum energy density, for large dimension gauge groups is cumbersome. Numerical preci- 
sion becomes a factor in the minimisation of the energy density. This problem is magnified 
in the calculation of tadpole improved results using the iterative procedure employed in 
Ref. 4. We thus abandon the use of variational wave functions here. Instead we make use 
of the one plaquette wave function of Eq. (3) and define a simple dependence of c(f3). Here 
we choose c((3) = (3. 

Many Hamiltonians are available for LGT calculations [22]. In Ref. 4 we used Kogut- 
Susskind [23], improved and tadpole improved Hamiltonians to calculate glueball masses for 
iV < 5. Here we employ the simplest of these, the Kogut-Susskind Hamiltonian, which is 
defined for pure SU(iV) gauge theory with coupling g 2 on a lattice with spacing a by 



To evaluate the matrix elements appearing in the glueball mass calculation we use the 




(4) 





where the plaquette operator is given by 




(6) 
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SU(iV) generating functions derived in Ref. 4, 

Gsu(N)(c,d) - 



dUe cTrU + dTrW 



SV(N) 



Z=— oo 



(7) 



and 



# m (c,d) = f dUexp [c(Tr/7 + Tr/7 f ) + dTr(/7 m )] 

JSV(N) 



E det 



Z=— oo 



fc=0 



(8) 



l<i,j<N 



Here the quantities inside the determinant are to be interpreted as the (i, j)-th entry of an 
N x N matrix. 

For the minimisation process, we use the same basis of rectangular states as in Ref. 4. 
With this basis the states J ++ and J , with J = or 2, are accessible. The J = 2 states 
require the use of a minimisation basis that contains no square states so that the excited 
states are invariant under rotations by n but not tt/2. The dependence of the massgaps on 
the coupling does not change significantly as iV is increased. We do however observe the 
appearance of what could possibly be additional low (3 scaling regions as iV is increased. 



B. Convergence with Z max 

In Ref. 4 the dependence of the variational parameter on the truncation of the /-sum 
in Eq. (7) was considered. The calculation of the variational parameter depends only on 
the plaquette expectation value. The dependence of the massgaps on the truncation should 
also be explored as their calculation incorporates additional matrix elements. We find that 
fast convergence is achieved as Z max is increased for the ++ massgaps. A typical example 
is shown in Fig. 1(a), where the massgaps are almost indistinguishable up to (3 = 80 on the 
scale of the plot. For the convergence is not as fast. A typical example is shown in 
Fig. 1(b). We find that in order to obtain convergence of the antisymmetric massgaps up to 
1/g 2 ?s 350 for SU(25) we need a truncation of Z max = 30. 
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(a) SU(4) AM++. 



(b) SU(3) AM—. 



FIG. 1: Example spin massgaps in units of g 2 /a demonstrating different convergence properties 
when truncating the £-sum of Eq. (7) at I = ±£ max . 

C. The small (3 minima 



The massgaps for iV > 2, in units of Ng 2 /a do not differ significantly as functions of 
1/g 2 . The lowest lying ++ glueball masses have the characteristic form shown in Fig. 1(a) 
for each N. We observe a minimum at small f3 and a scaling plateau at large (3. 

In this section we examine the minima occurring in the lowest ++ and glueball 
masses and consider the possibility that they may correspond to a scaling region. It appears 
that for SU(3) this region has been interpreted as a scaling region by some authors [24- 
26]. From the scaling arguments of Ref. 4, by changing the functional dependence of the 
parameter appearing in the vacuum state the scaling region can be modified. In this way 
it is possible to choose the variational parameter's dependence on [3 such that the region 
in which the mass gap takes its minimum is stretched over a large j3 interval. In the same 
way, the scaling plateaux occurring for large (3 may be compressed. With this in mind, it 
is unclear which scaling region we should take to be the correct one. For this reason we 
analyse both potential scaling regions. We start, in this section, with the small j3 minima 
and consider the large (3 plateaux in the next section. 

Fig. 2 shows the lowest lying SU(iV) ++ glueball masses in 2+1 dimensions in units of 
g 2 N/a as functions of 1/g 2 for various N e [3,25]. It is apparent that the minima depend 
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only weakly on N. We see that the minima appear to approach a finite limit from below 
as iV — > oo. On the scale of Fig. 2 the minima corresponding to N = 15 and N = 25 are 
barely distinguishable. 

Let us denote the minima in Fig. 2 by aAM c / (Ng 2 ) and consider them as a functions of 
N. Fitting aAM c /(Ng 2 ) to the model 

^ = 7l + ^_ (9) 

N g2 71 + ^ W 

for > 8 gives best fit parameters 

71 = 0.83262 ±0.00022 

72 = -0.97 ±0.13 

73 = 0.990 ±0.035. (10) 

We fit on the data iV > 8 to minimise pollution from next to leading order corrections. If we 
assert that the power of the leading order correction for finite iV must be an integer power 
of 1/N, then from Eq. (10) that order must be 2 in agreement with the predictions of large 
N QCD. Fitting to a model with 1/N 2 corrections for N > 8 gives 

W (^) = 0.83256±0.00007- a9753 ^ a ° 072 . (11) 

We can go further and attempt to find the next to leading order corrections by fitting to 
the model 

^ = 7l + ^ + ^. (12) 
Ng 2 N 2 N 2 ^ 1 ; 

The best fit parameters when fit on the whole data set are 

71 = 0.83287 ±0.00011 (13) 

72 = -1.116 ±0.028 (14) 

73 = 1.715 ±0.085 (15) 

74 = 1.627 ±0.050 (16) 

which is in disagreement with integer power next to leading order corrections. A good fit 
however is achieved to the model with 1/N 4 next to leading order corrections: 

0.986 ±0.010 2.772 ±0.090 
fi 2 (N 2 ) = 0.83233 ± 0.00018 — + ^ . (17) 



The precise locations of the minima, aAM c / (Ng 2 ), are plotted as a funtion of 1/N 2 in Fig. 4 
along with the fitted models of Eqs. (11) and (17). 

From Eq. (11), our best fit result for the N -> oo limit of aAM c /(Ng 2 ) is 0.83256 ± 
0.00007. The error here is purely statistical. We should, as always, expect a significant 
systematic error attributable to our choice of ground state and minimisation basis. This 
result may be compared to the Monte Carlo calculation of Lucini and Teper who obtain the 
N — > oo limit of the lowest ++ glueball mass as 0.8116 ± 0.0036 [2]. This is comparable to 
the result presented here. The result of Lucini and Teper is based on a linear extrapolation 
to the N — > oo limit of 2 < N < 6 data. Their result differs significantly from the one 
presented here at the leading order finite N corrections. The correction term obtained by 
Lucini and Teper, —(0.090 ± 0.028)/iV 2 , has the same sign but is significantly smaller than 
the one presented here. When our data is fit on the range 3 < iV < 6 (we don't obtain a small 
f3 minimum for SU(2)) our slope is halved but is still significantly larger than that of Lucini 
and Teper. It should be pointed out that Lucini and Teper's calculation was performed in 
the Lagrangian approach in which the coupling is the so called Euclidean coupling qe- The 
Euclidean coupling and the Hamiltonian coupling, g 2 , are equal up to order g 2 E corrections. 
The precise relation between the couplings is only known for small g 2 E [27], a case which 
does not apply for the small (3 minima. It appears that our simple calculation induces a 
level crossing whereby the lowest mass glueball state switches to a higher mass state beyond 
the small (3 minimum. Although we do not have an explanation for this, presumably, by 
including additional states in the minimisation basis or implementing a more complicated 
vacuum state, the level crossing would no longer appear and the small (3 minima would 
extend into large (3 scaling regions. This should be checked in further studies. A glueball 
mass extracted from a large (3 scaling region can be confidently compared to a corresponding 
Lagrangian calculation, for in that region of couplings the ratio of g 2 E to g 2 is unity. 

We obtain small (3 minima for the massgaps but only for N > 5. Plots of these 
massgaps with N > 5 are shown in Fig. 3. The behaviour is significantly different to 
that observed for the symmetric massgap. By N = 25 the minima do not appear close to 
convergence. Indeed if convergence is occurring at all, it is significantly slower than was 
observed for the ++ state. To explore this further, in Fig. 5 we plot the minima of Fig. 3, 
which we again denote by aAM c / (Ng 2 ), as a function oil/y/N. We observe linear behaviour 
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in the large N limit. Fitting the N > 9 data to the model 



aAM c 72 , . 

TvF = 71 + iv^' (18) 



gives best fit parameters 



71 = 0.41 ±0.01 

72 = 1.25 ±0.02 

7s = 0.98 ±0.04, (19) 

which is consistent with 73 = 1. Fitting the N > 9 data to a model with leading order finite 
iV corrections starting at 1/ \/N gives 

MN) = 0.41390 ± 0.00007 + L255 ^_°- 003 , ( 20 ) 



As was done for the symmetric case we can attempt to go further and find the form of the 
next to leading order finite N corrections. Fitting the N > 5 data to the model 

, 2W = 7l + ^ + _*_. (21) 



gives best fit parameters 

71 = 0.410 ±0.001 (22) 

72 = 1.2718 ±0.0045 (23) 

73 = -76 ±23 (24) 

74 = 9.5 ±0.4 (25) 

indicating vanishingly small next to leading order corrections for iV > 5. The models v\ and 
u 2 are plotted against the N > 5 data in Fig. 5. 

The small j3 minima appear to converge to a non-zero N — > 00 limit. The convergence is 
significantly slower than for the small (3 minima in the symmetric case. It is also clear that 
the leading order finite N corrections are not the expected 0(1/N 2 ) but rather 0(l/y/N). 
The N — >■ 00 limit of 0.41 ± 0.01 is approximately half the corresponding ++ result. Such 
a state does not appear in the calculation of Teper. 

It is possible that the small (5 minima are spurious. Dimensional analysis gives the 
expected scaling form but does not allow us to decide which scaling region is preferable, or 
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FIG. 2: L max = 4 SU(iV) lowest mass symmetric massgaps in 2+1 dimensions in units of g 2 N/a 
as functions of l/g 2 . The l-sum of Eqs. (7) and (8) truncated at l max = 5. 

if one is a lattice artifact. However, based on lattice calculations to date, we expect the ++ 
glueball to be lighter than the in the continuum limit. Spurious scaling regions have 
been observed in other lattice calculations [28]. It is clear that additional analysis is needed 
here. Again, an important step will be to include additional non-rectangular states in the 
minimisation basis. 

D. The large (3 plateaux 

Having considered the small j3 minima as possible scaling regions we now move on to 
the large /3 plateaux. This scaling region appears for all N and for all states considered 
and so its interpretation as a genuine scaling region is less dubious. Furthermore, the 
glueball mass results extracted from these scaling regions may be confidently compared to 
the corresponding Lagrangian calculations since the ratio of the Hamiltonian to Euclidian 
coupling is unity up to small 0(g 2 E ) corrections. 

We start with the ++ state for which the best scaling behaviour is obtained. We calculate 
the five lowest lying massgaps corresponding to the five lowest glueball masses accessible 
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FIG. 3: L max = 6 SU(iV) lowest mass anti-symmetric massgaps in 2+1 dimensions in units of 
g 2 N/a as functions of l/g 2 . The l-sum of Eqs. (7) and (8) truncated at l mSLX = 5. 
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FIG. 4: The continuum limit lowest mass symmetric 2+1 dimensional massgaps in units of Ng 2 /a 
as a function of 1/N 2 taken from the small j3 minima. The dashed lines are fits to Eqs. (11) and 
(17). 
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FIG. 5: The continuum limit lowest mass antisymmetric 2+1 dimensional massgaps in units of 
Ng 2 /a as a function of \/y/~N taken from the small (3 minima. The dashed lines are fits to 
Eqs. (20) and (21). 

with our choice of ground state and minimisation basis. We consider values of N in the 
range 3 < iV < 25. For each N considered we find a large j3 scaling plateau for each of 
the five lowest mass states. The lowest lying massgap is shown in Fig. 6 for a range of N. 
We observe that in units of Ng 2 /a the massgaps do not depend strongly on N and that in 
the scaling region they appear to approach a finite limit. Similar observations can be made 
for the four higher mass states obtained although the scalaing behaviour is less precise. We 
show the second lowest glueball mass in Fig. 7 as an example. A continuum limit for each 
massgap is obtained in the scaling region by fitting to a constant. For each fit we use a 
region of at least 10 data points which minimises the standard error. The continuum limit 
results for the five lowest lying massgaps, denoted by aAM c /(Ng 2 ), are shown as functions 
of 1/N 2 in Figs. 6 and 9. Also shown in the plots are fits to models with leading order large 
iV corrections starting 1/iV 2 , 



++ — _i_ — _i_ Ti 

K i ~ Pi + N 2 + N 4 ■ 



(26) 



Here the subscript i labels the i-ih lowest glueball mass. The values of the best fit parameters 
are given in Table I. 
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TABLE I: The best fit parameters for the five lowest energy ++ massgaps when fitting Eq. (26) 
to the available data. 



Pi 



1 

2 
3 
4 
5 



1.23526+0.00023 
2.36924±0.00072 
2.88446±0.00071 
3.35422±0. 00047 
3.7667±0.0013 



-1.540±0.018 
-2.478±0.047 
-3.435±0.017 
-3.476±0.012 
-4.114±0.092 



1.97±0.16 
-1.62±0.42 

Q a 
Q a 

1.88±0.83 



"Set to zero to obtain a stable fit. 



These results should be compared with those of Lucini and Teper who obtain large N 
glueball masses in units of Ng 2 /a in the ++ sector with linear fits given by 

0.090(28) 



++ 

++* 

Q + + ** 



0.8116(36) - 
1.227(9) - 
1.65(4) - 



N 2 

0.343(82) 

N 2 
2.2(7) 



N 2 



(27) 



The fit for the ++ ** state is obtained using 4 < iV < 6 data. The remaining fits are obtained 
using 2 < iV < 6 data. We find that the lowest glueball mass extracted from our large (3 
plateaux is consistent, in the N — ■> oo limit, with the state which Lucini and Teper label 
++ **. The slopes of the fits are of the same sign but differ significantly in magnitude. The 
q++** gt a ^ e f Lucini and Teper does not appear in our data in the form of a large /3 scaling 
plateau. There is however a hint of an approach to scaling in the vicinity of their result 
in our second lowest massgap data. This effect is only visible in our data for iV > 13 and 
occurs for small (3 as shown in Fig. 10. The values of (3 for which this effect appears are quite 
close to those for which the small beta minima are observed in the lowest mass eigenstate. 
Similar effects are observed in the data but the effect is much less convincing with the 
data currently available. 

Let us consider this small (3 region as a possible scaling region. We fit a constant to the 
available N > 15 data on a range of at least 4 data points chosen so that the standard error 
is minimised. These scaling values are plotted as a function of 1/N 2 in Fig. 11. Also shown 
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is the best fit linear model 

5.83 ±0.97 , x 

K ++ = 1.7605 ±0.0032 — . (28) 

This produces an A" — > oo limit which is close to but inconsistent with the result obtained 
by Lucini and Teper for their ++ ** state. 

We now move on to the states. For these states the scaling behaviour is not as 
precise as that obtained for the ++ states. Again we calculate the five lowest lying massgaps 
corresponding to the five lowest mass glueballs accessible with our choice of minimisation 
basis and ground state. For each A" considered up to 25 we find a large (3 scaling plateau 
for each of the five states. The lowest lying massgap is shown in Fig. 12 for a range of N. 
The second lowest energy massgap is shown in Fig. 13. Considerably less data has been 
obtained for the states due to the large Z max required for convergence as discussed in 
Section III B. Despite this, the results qualitatively replicate those of the ++ states. We 
observe that in units of Ng 2 /a the massgaps do not depend strongly on A" and that in the 
scaling region they appear to approach a finite limit. Similar observations can be made for 
the four higher mass states although the scaling behaviour is less precise. The continuum 
limit values extracted are plotted as functions of 1/N 2 in Figs. 14 and 15. The dashed lines 
show fits of the continuum limit values to the model 

+ + (29) 

with the parameters for each excited state given in Table II. 

These results should be compared again with those of Lucini and Teper who obtain large 
A" glueball masses in units of Ng 2 /a in the sector with linear fits given by 

0.14(20) 



0"" : 1.176(14) + 
0""* : 1.535(28) 



N 2 
0.35(35) 

A^ 2 



0—: 1.77(13) ± °" 24 A ( i 61) - (30) 

N 2 

The result for ** is from Ref. 1. All fits were obtained using 3 < A^ < 6 data. The 
A^ — > oo limit of our lowest lying massgap extracted from the large f3 plateaux is consistent 
with Teper's ** result. 

Having considered spin states we now move on to spin 2, the only other spin accessible 
when using a basis of rectangular states. For the case of spin 2 the scaling behaviour of 
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TABLE II: The best fit parameters for the five lowest energy massgaps when fitting Eq. (29) 
to the available data. 



Pi 



(It 



1 

2 
3 
4 
5 



1.8896+0.0011 
2.8930±0.0044 
3.20871±0.0047 
3.83647±0.0024 
4.0759±0.0019 



-1.829±0.068 
-2.96±0.23 
-3.35±0.22 
-3.41±0.05 
-2.45±0.12 



6.95±0.62 
6.6±2.0 
8.8±1.9 

Q a 

-7.4±1.1 



"Set to zero to obtain a stable fit. 



the massgaps is significantly worse than for spin 0. More troublesome is the fact that 
the convergence of the 2 ++ massgaps with increasing Z max is slower than the case of . 
The situation for the 2 state is markedly better with the convergence of massgaps with 
increasing / max being no different to that of ++ . For this reason we concentrate on the 2~~ 
sector here. 

The 2 states produce less precise scaling behaviour than the spin states. However 
large (5 plateaux appear for each of the five lowest lying massgaps. We use these regions to 
estimate their respective continuum limits. For these states small f3 minima do not appear. 
The lowest lying 2 massgap is shown in Fig. 16 for a range of N. Once again we observe 
that in units of Ng 2 /a the massgaps do not depend strongly on iV and that in the scaling 
region they appear to approach a finite limit. Similar observations can be made for the four 
higher mass states obtained, although the scaling behaviour worsens as the mass of the state 
increases. The estimated continuum limit values extracted are plotted as functions of 1/N 2 
in Figs. 17 and 18. The dashed lines show fits to the model 



a— _ 3i_ Jj_ 
* Pi + N 2 N 4 ' 



(31) 



The best fit parameters for each excited state are given in Table III. 

Once again these results should be compared with those of Lucini and Teper who obtain 
large N glueball masses in units of Ng 2 /a in the 2 sector with linear fits derived from 
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TABLE III: The best fit parameters for the five lowest energy 2 massgaps when fitting Eq. (31) 
to the available data. 



i Pi Qi n 

1 3.20379±0.00006 -2.9571±0.0032 5.73912±0.02585 

2 4.07376±0.00066 -2.189±0.050 -9.46±0.43 

3 4.96259±0.00090 -4.947±0.045 1.99±0.35 

4 5.26717±0.00040 -6.221±0.031 15.92±0.26 

5 5.744±0.024 -5.13±0.56 a 



a Set to zero to obtain a stable fit. 
3 < N < 6 data given by 

9" ,fii^ °- 1Q ( 42 ) 
2 : 1.615(33) j^- 

2— 1.87(12) (32) 

The 2 * result is from Ref. 1. The N — > oo results presented here are significantly higher. 
This time no correspondence between our results and those of Lucini and Teper can be 
obtained. 

We finish this section by presenting the mass spectra obtained for the ++ , and 2~~ 
sectors. The results are summarised by the plots in Figs. 19 — 21. 

E. An empirical observation 

In this subsection we present an empirical observation that could possibly be useful in 
the construction of simple glueball models. Since Lucini and Teper include many small area 
Wilson loops in their construction of excited states, we assume that they calculate the lowest 
few glueball states without omission. We assume the calculations presented here give higher 
mass glueball states. Presumably, the correct enumeration of these states is not possible 
without the inclusion of nonrectangular states in the minimisation basis. For example, the 
third lowest mass state in a given sector calculated here could be the true seventh lowest 
mass state. However, when our results are combined with those of Lucini and Teper an 
interesting empirical observation can be made. We can choose a labelling of the ++ , 0~~ 
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FIG. 6: Lowest mass symmetric 2+1 dimensional massgaps in units of Ng 2 /a as a function of l/g 2 
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FIG. 8: Continuum limit lowest mass ++ SU(iV) massgap in units of Ng 2 /a as a function of 1/N 2 . 
The dashed line is the fit to the quadratic model of Eq. (26). 

and 2 excited states such that the large N limit of their masses lie on a straight line as 
shown in Fig. 22. We have included the results of Lucini and Teper for , *, 2 and 
2 * which do not correspond to any of the states calculated here. The straight lines are 
fits to the model, 



m n (J PC ) = 7l (2n + 72 ) 



(33) 



where 71 and 72 are parameters and 72 is restricted to integer values. For the J PC states 
considered we obtain the following best fit models 

m n (0++) = (0.256 ± 0.002)(2n + 1) 
m„(0"") = (0.151±0.002)(2n + 5) 

m n (2— ) = (0.1495 ±0.0008)(2n + 7). (34) 

Here n > 1 labels the n-th lowest mass state. We note that the fit improves in accuracy 
as n is increased. It is interesting to note that the slopes of the and 2 models are 
consistent suggesting that the constant of proportionality does not depend on J. Another 
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FIG. 9: Continuum limit ++ SU(iV) massgaps in units of Ng 2 /a as functions of 1/N 2 . The dashed 
lines are fits given in Eq. (26). 

interesting observation can be made by recasting the models in the form 

m n (0++) = (0.256 ±0.002)(2n + + 1) 
m n (0"") = (0.151 ± 0.002)(2n + 4 + 1) 

m„(2~) = (0.1495 ±0.0008)(2n + 6 + 1). (35) 
We notice a similarity with the two dimensional harmonic oscillator spectrum, 



E n oc 2n + J + 1, 



(36) 
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FIG. 14: The 2+1 dimensional continuum limit lowest SU(iV) glueball mass in units of Ng 2 /a 
as a function of 1/N 2 . The dashed line is the fit to the quadratic model of Eq. (29). 

when we take into consideration the ambiguity modulo 4 of spin identification on the lattice. 
With this in mind we propose a simple model for the J PC spectrum 

m n {J PC ) = Cpc [2n + i{J PC ) + l] , (37) 

where Cpc is a spin independent parameter and 7(J PC ) is an integer for which 7(J PC ) = 

J mod 4. From Eq. (34) we have C ~ 0.15 and C++ = 0.256 ± 0.002. To check this simple 

model we can attempt to predict the lowest lying states obtained by Lucini and Teper in 
the sectors that have not been considered in this paper. We start with the 2 ++ sector in 
which Teper and Lucini obtain 1.359(12) and 1.822(62) for the N -> oo limit of the 2++ 
and 2 ++ * glueball masses in units of Ng 2 /a. The predictions of Eq. (37), with 7(2 ++ ) = 2 
and C++ = 0.256, are 1.28 and 1.79 for the 1st and 2nd lowest glueball masses in units of 
Ng 2 1 a. As we would expect the prediction is better for the higher mass state. We can also 
consider the 1 ++ sector which has an equivalent spectrum to 1~ + due to the phenomenon 
of parity doubling [1]. Lucini and Teper obtain 1.98(8) for the mass of the 1 ++ glueball in 
units of Ng 2 /a. The model of Eq. (37), with 7(1++) = 5 and C++ = 0.256 gives 2.048. In 
the 1 sector the agreement is not as good, with Lucini and Teper obtaining 1.85(13) for 
the mass of the 1 glueball and the model of Eq. (37), with 7(1 ) = 9 and C = 0.15, 
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FIG. 15: Continuum limit SU(iV) massgaps in units of Ng 2 /a as functions of l/N 2 . The 
dashed lines are fits given in Eq. (29). 



giving 1.812. 

It will be interesting to see if the model suggested here stands up to further calculations 
with an extended minimisation basis and in other J PC sectors. The analytic techniques pre- 
sented here are at a great advantage to the standard Monte Carlo techniques of Lagrangian 
LGT for the purpose of testing glueball models. High order excited states are readily acces- 
sible; in a variational calculation with s states in the minimisation basis, s glueball states 
are accessible. This is in contrast to competing Lagrangian calculations in which only 3 
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as a function of \/N 2 . The dashed line is the fit to the quadratic model of Eq. (31). 
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FIG. 18: The 2+1 dimensional continuum limit 2 SU(iV) massgaps in units of Ng 2 /aas functions 
of 1/N 2 . The dashed lines are fits given in Eq. (31). 



states are currently accessible in some J PC sectors. 



F. Discussion 

In this section we have calculated variational mass spectra for pure SU(iV) gauge theory 
with a simple basis of rectangular states in 2+1 dimensions. Such a basis is easy to work 
with computationally and is therefore a good starting point. However to accurately explore 
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the pure gauge mass spectra additional states need to be included. Perhaps a more suitable 
basis would be the set of all states with an area less than some maximum value which would 
define the order of the calculation. The next stage beyond this would be to include higher 
representation states in the form of Wilson loops in which some links are covered more than 
once. In this way multiple glueball states could presumably be explored. The inclusion of 
additional small area states in the basis, in particular states which are not symmetric about 
reflections in coordinate axes, would also allow the calculation of spin 1 glueball masses and 
also massgaps in the P = —C sector. 

It is important to point out that the inclusion of additional states in our minimisation 
basis does not present a significant challenge. The only complication would be in counting 
the possible overlaps of particular states. While more difficult than the counting required 
here, the process could presumably be automated using symbolic programing and techniques 
from graph theory. The analytic techniques used here would still be applicable until higher 
representation states were included in the minimisation basis. At that stage further character 
integrals would be required to handle the additional integrals arising in the calculation. 

A final extension of the method presented here would be to make use of an improved 
ground state; a ground state which includes additional Wilson loops in the exponent. How- 
ever, the advantage of using analytic techniques would then be lost, unless new technology 
was developed for tackling the resulting integrals. One would need to calculate the re- 
quired matrix elements via Monte Carlo simulation on small lattices, losing the advantage 
of working on an infinite lattice with analytic expressions that we have here. 

IV. CONCLUSION 

In this paper we have applied the analytic techniques developed in Ref. 4, in a study of 
the large N glueball mass spectrum in 2+1 dimensions. 

In Section III we calculated glueball masses at finite N, with N as large as 25 in some 
cases. This allowed accurate N — > oo extrapolations to be made. Evidence of leading order 
1/N 2 finite N corrections to the glueball masses was obtained, confirming a specific predic- 
tion of large N gauge theory. The interpretation of the possible scaling regions was discussed 
and agreement with the Lagrangian study of Lucini and Teper was obtained in some cases 
in the N — > oo limit. The discrepancies are attributable to our use of only rectangular states 
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in the minimisation basis. Without the inclusion of small area, nonrectangular states, it is 
likely that some low energy states are inaccessible. Further work will fill in the incomplete 
spectra presented here. Interesting empirical observations were made when the results pre- 
sented here were combined with those of Lucini and Teper; the enumeration of excited states 
can be chosen so that the glueball mass spectrum has the structure of a two dimensional 
harmonic oscillator. To develop this observation into a model would require a more complete 
calculation of the mass spectrum with nonrectangular states in the minimisation basis and 
possibly a more complicated vacuum trial state. 

Clearly much work remains to be done before an accuruate picture of the 2+1 dimensional 
pure SU(iV) gauge theory spectrum is achieved within a Hamiltonian variational approach. 
We have demonstrated however that the analytic techniques of Hamiltonian LGT can be used 
for massgap calculations with iV as large as 25 on a desktop computer. This is significantly 
closer to the N — > oo limit than is currently possible with the Lagrangian approach using 
supercomputers. The inclusion of additional states in the minimisation basis, while not 
presenting a major challenge, will allow a thorough study of the 2+1 dimensional pure 
SU(iV) mass spectrum at least up to N = 25. 
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